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Abstract 

We present experimental results for spinodal decomposition in polymer mixtures of gelatin and 
dextran. The domain growth law is found to be consistent with t^/^-growth over extended time- 
regimes. Similar results are obtained from lattice simulations of a polymer mixture. This slow 
growth arises due to the suppression of the bulk mobility of polymers. In that case, spinodal 
decomposition is driven by the diffusive transport of material along domain interfaces, which gives 
rise to a t^/^-growth law. 
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A homogeneous binary mixture (AB) becomes thermodynamically unstable when it is 
quenched into the miscibihty gap. The subsequent evolution of the mixture is characterized 
by the emergence and growth of domains enriched in either component This 
process of domain growth or coarsening or spinodal decomposition is of great technological 
and scientific importance, and is relevant in diverse applications, e.g., the food-processing 
industry, metallurgy, materials science, structure formation in the early universe, etc. The 
standard experimental tools to characterize domain growth are (a) the time-dependence 
of the typical domain size and (b) the functional form of the structure factor, which 
depends on the system morphology. The domain growth law depends critically on the 
mechanism which drives phase separation. For example, diffusive transport of material 
through bulk domains yields the Lifshitz-Slyozov (LS) growth law ~ ^^^t^i 'whereas 
diffusion along surfaces (interfaces) gives rise to a slower growth law ^{t) ~ t^/^ 
For binary fluid mixtures, the intermediate and late stages of phase separation are driven by 
the advective transport of material. The relevant growth laws (in (i = 3) are £{t) ~ t in the 
viscous hydrodynamic regime, followed by £(t) ~ t^/^ in the mertial hydrodynamic regime 

HQS. 

In the food-processing industry, many applications involve domain growth in demixing 
polymer blends or quasi-binary polymer solutions (i.e., two kinds of polymers in a solvent). 
Polymeric phase separation has been studied experimentally in films of polymer blends 8|, 
y. nil III, in thicker samples Q, and in soluted mixtures The typical time-scales 
of decomposition in these experiments vary from seconds Q| to days To characterize 
the speed of the coarsening process, one often fits a power-law to the location gmax of the 
maximum in the structure factor, as a function of time. Measurements of the resulting 
growth exponent a in g^ax ~ range from to 0.45 0, 10, ^|. Wilzius and Gumming 
12 1 report a value of a ~ 0.28; Kuwahara et al report an evolution of this exponent 



from a ~ to a ~ 0.3; and Takeno and Hashimoto ^ measure a value between 0.25 and 
0.33. Clearly, there is no consensus on the experimental side regarding the growth law. 
Moreover, the dominant mechanism for domain growth may differ in various experiments 
and time-regimes, explaining the diverse exponents reported in the literature. 

In this letter, we present results from an experimental and numerical study of spinodal 
decomposition in polymer mixtures in the high- viscosity regime (i.e., without hydrodynam- 
ics). First, we will present results from a light-scattering study of phase separation in a 
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mixture of gelatin and dextran, soluted in water. For deep quenches, domain growth in 
this system is consistent with the growth law i(t) ~ t^^^. Next, we will present numerical 
results obtained from a lattice simulation of spinodal decomposition in a polymer mixture. 
Again, the domain growth law is found to be compatible with i(t) ~ t^^^ over an extended 
time-regime. We interpret these growth laws in the context of phase separation driven by 
surface diffusion, as the bulk mobility is drastically reduced in these polymer mixtures due 
to a large free-energy barrier for releasing a component from its concentrated phase. We 
conclude this letter with a discussion of the growth exponent which characterizes surface 
diffusion, and the crossover from surface to bulk diffusion. 

Let us first present details of the experiments. The chemicals used for the spinodal 
decomposition experiments were gelatin (fish gelatin with a high molecular weight, Multi 
Products, Lot 9187, box 47) and dextran (obtained from Sigma, industrial grade with an 
average molecular weight of 181 kg/mol, D-4876, Lot 41K1243). These chemicals were used 
without further purification. The concentration in the final sample for gelatin was 2.46% 
and for dextran 2.40% w/w. No salt was added to the solution. However, azide was added 
to the stock solutions to halt bacterial growth. The light-scattering set-up used is shown in 
Fig. 2 of Tromp et al A He/Ne laser (633 nm) was used to generate a scattering pattern 
that was projected onto a screen of tracing paper. This is recorded by a CCD-camera (EG 
& G PARC, Model 1430P). (See also Lundell et al |l6|.) The usable g-range was 0.3-7 fim-\ 

The gelatin-dextran mixture described above is homogeneous at room temperature and 
starts to demix at a temperature of ~ 8°C. The phase separation is clearly visible at 6°C. 
This is close to the gelation temperature of the gelatin, which is estimated as 8-10°C Q|. 
We specifically selected this mixture to obtain a high viscosity, so that hydrodynamic effects 
can be ignored on the time-scales of our experiments. 

Starting with a homogeneous mixture, we performed either a shallow quench to a final 
temperature of 5°C, or a deep quench to 1°C. From then on, the structure of the coarsening 
mixture is recorded at regular times with light scattering. Specifically, we obtain the typical 
domain size from £(t) = 2n / q^g^^it) , where qmaxit) is the location where the spherically- 
averaged structure factor S{q,t) reaches its maximum. Due to practical reasons, the time 
to of the onset of phase separation is not known precisely in our experiments. We therefore 
determine to from our measurements of S{q, t) as follows. Assuming that the domain growth 
process is described by gmax(t) ~ (t — to)~" with a ~ 1/4 or 1/3, we perform a linear least- 
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squares fit in which gmax is plotted vs. t; the time of zero-crossing is then our estimate of 
to- Since it takes some time before the structure factor develops a clear peak, we repeated 
this procedure also for the locations q'^^^ and q'^^^ of the maxima in qS{q,t) and q^S{q,t). 
This results in six estimates for Iq, which we average. We use the average value of to to 
make a double- logarithmic plot of ^{t) vs. (t — to); the result is shown in Fig. ^ 




FIG. 1: Time-dependence of length scale from experiments, plotted on a log-log scale. We study 
a deep quench to 1°C (left frame) and a shallow quench to 5°C (right frame). The lines denote the 
best linear fits to the experimental data. The corresponding growth exponents are specified in the 
figure. 

As the straight lines in the figure indicate, domain growth after a deep quench to 1°C 
is fitted well by £(t) ~ (t - to)°■^^ and after a shallow quench to 5°C by £(t) - (t - tof-'^'^. 
The deep quench shows convincingly a domain growth exponent in agreement with surface 
diffusion. The shallow-quench experiment is not so conclusive. The larger growth exponent 
indicates that both surface and bulk diffusion contribute to coarsening, and we are in a 
crossover regime between a = 1/4 and a = 1/3. As we will discuss shortly, the bulk 
mobility is enhanced at higher temperatures where A-rich bulk domains contain significant 
amounts of B, and vice versa. 

We have also performed extensive simulations of spinodal decomposition in a quasi-binary 
polymer mixture. The specific model that we used is a lattice polymer model. Polymers with 
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bonds are described in this model as self-avoiding and mutually-avoiding walks on a face- 
centered-cubic (FCC) lattice. The polymers move by a sequence of single-monomer jumps to 
nearest-neighbor lattice sites. For monomers which are neighbors along the same polymer, 
the steric interactions are lifted (and a site can thus be occupied by two or more adjacent 
monomers), as this allows enhanced reptation. Both reptation and Rouse dynamics are 
captured qualitatively in this model. Therefore, it can be used to study dynamical properties 
with some degree of realism, as long as hydrodynamic interactions are not essential. The 
model is described in detail by Heukelum and Barkema It has previously been used to 
study the fractionation of polydispersed polymer mixtures 13, diffusion and exchange of 
adsorbed polymers [20], and desorption of polymers 

The current study involves 128,000 polymers, each consisting of = 50 bonds and 
connecting 51 monomers. These are located on a 3-dimensional FCC lattice with 8 x 10^ 
sites and periodic boundary conditions. Half of the polymers are labeled as type A, and the 
other half as type B. Polymers of different type repel each other: the total energy is equal 
to the number of pairs of neighboring sites occupied by different-type polymers, multiplied 
by an energy scale J. If an attempted Monte Carlo (MC) move increases the total energy 
by /\E, the acceptance probability of this move is Pa = exp{—j3AE), where (3 = l/lksT) is 
the inverse temperature. Moves which do not raise the total energy are always accepted. 

The simulations start from a mixed state, generated as an equilibrium configuration at 
high temperature (/5 = 0). At t = 0, the temperature is quenched to either pj = 0.05 or 
/3J = 0.1. At these low temperatures, the repulsive AB interactions induce phase separation 
into domains rich in either A or B. The size of these domains grows in time - evolution 
snapshots for a typical run are shown in Fig. |21 To characterize domain growth, we assign 
to lattice site i occupied by an A-type (B-type) polymer the value cxj = -|-1 (a, = —1). The 
domain size is defined as the distance £ of the first zero-crossing of the spatial correlation 
function of (Xj, and we measure it at various times. The time-dependence of i{t) is shown in 
Fig. El After an initial transient regime, both data sets (for (3 J = 0.05,0.1) are consistent 
with the growth law i{t) ~ t^/^ over extended time-regimes, viz., ~ 2 decades for pj = 0.05 
(denoted by circles), and ~ 2.5 decades for pj = 0.1 (denoted by squares). 

The above experimental and numerical results can be understood in the following theo- 
retical framework. The kinetics of phase separation is described by the Cahn-Hilliard (CH) 
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FIG. 2: Evolution pictures from simulations of spinodal decomposition in polymer mixtures. 
Details of the simulation are described in the text. The quench temperature was f3J = 0.1, and 
the system size is 200^ The species A is marked in black, and the species B is not marked. The 
time in Monte Carlo steps (MCS) is indicated above each snapshot. 
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FIG. 3: Time-dependence of length scale from simulations, plotted on a log-log scale. The time 
is measured in MCS. We show data for quenches to (53 = 0.05 (circles), and /3J = 0.1 (squares). 
The lines correspond to the growth law £(t) ~ 



model 
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dt 
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where ip{r,t) is the order parameter at space point r and time t. Typically, ip{r,t) ~ 
PA{r,t) — pB{r,t), where pa and pb denote the local densities of species A and B. We have 
neglected thermal fluctuations in Eq. ((H), and considered the general case of a ^-dependent 



mobility D{ijj) 22; |23|, |2J| . Further, F[iIj] denotes the Helmholtz potential, which is usually 



taken to have the ip^-foim: F[ip] c::^ j df —iIj^/2 + iIj'^/4+ {Vtp)'^/2 , where we have used 
the order-parameter scale and the bulk correlation length to express F[ip] in dimensionless 
units. 

Let us consider a situation where the mobility at the interfaces {ip = 0) is Dg, and that in 
the bulk {ip = ±1) is Dh, with Dh < Ds- This difference in interfacial and bulk mobilities can 
result from low-temperature dynamics, where there is a high energy barrier for an A-atom 
to enter a B-rich bulk domain and vice versa. Alternatively, physical processes like glass 
formation j^] or gelation [3, can also result in a drastic reduction in bulk mobility. 
A simple functional form for the mobility in the above case is D{iIj) = Ds{l — aip"^), a = 
1 — Dfj/Dg. Notice that < a < 1 for Df, < Dg. Thus, Eq. can be written as 



(1 - ai)^) V {-^ + - V^^) 



(2) 



where we have absorbed Dg into the time-scale. 
d 



dt 



ijj{f, t) = {l- a)W\-^ + - V^^) + aV • (1 - V {-^ + - V^^) 



(3) 



where the first term on the RHS corresponds to bulk diffusion. This term is absent in the 
limit a = 1 {Db = 0). The second term on the RHS corresponds to surface diffusion. It is 
only relevant at interfaces, where ~ 0. The location of the interfaces rj(t) is defined by the 
zeros of the order-parameter field, ip [fi{t),t] = 0. Let us focus on a particular interface in the 
(i-dimensional case. We denote the normal coordinate as n and the interfacial coordinates 
as a [with dimensionality {d — 1)]. Then, the normal velocity Vn{3,t) obeys the integro- 
differential equation 



aality 

m 



4 j da'G[fi{d),fi{a')]vnia',t) ^ (1 - a)aK{d,t)+Aa j da'G[fi{d),fi{a')]V^K{a',t). (4) 

Here, K{a, t) is the local curvature at point a on the interface, and a is the surface tension. 
The Green's function G{x^x') is obtained from — V^G(x, a;') = 5{x — x'). 

A dimensional analysis of Eq. (jU yields the growth laws due to surface and bulk diffusion. 
We identify the scales of various quantities in Eq. Q as 



[da\ 



d£ 
di' 



[G] 



02-d 



(5) 



This yields the crossover behavior of the length scale as 

i{t) ~ t<^t„ 

~ [(1 - i>ie, (6) 

where tc ~ ^^/[(l — Notice that the asymptotic regime obeys the LS growth law, 

which is characteristic of growth driven by bulk diffusion. However, the crossover to the 
LS regime can be strongly delayed if the bulk mobility is suppressed relative to the surface 
mobility, i.e., D^/ Dg ^ or a ^ 1. 

In conclusion, we have presented results from both experiments and simulations of spin- 
odal decomposition in polymer mixtures. In both cases, we find an extended time-regime 
where the domain growth law is compatible with £{t) ~ t^^^. This slower growth law arises 
because there is a strong suppression of mobility in bulk domains due to a large energy 
barrier to dissolve in the polymer-poor phase. Therefore, domain growth is primarily 
driven by diffusive transport along interfacial regions, which gives rise to a i^/^-growth 
law. Understanding the mechanism of demixing is not only important for identifying the 
domain growth exponent, it is also a pre- requisite for obtaining control over the structure 
formation. We hope that the results presented here will stimulate further experimental and 
numerical interest in this fascinating problem. 
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